INTERNAL AND EXTERNAL AXIAL CORNER FLOWS 

By Paul Kutler, 

NASA Ames Research Center 

Vijaya Shankar, Dale A. Anderson, 
Iowa State University 

• and Reese L, Sorenson 
NASA Ames Research Center 


SUMMARY 


The inviscid, internal and external axial corner flows generated by two 
Intersecting wedges traveling supersonically are obtained by use of a second- 
order shock-capturing, finite-difference approach. The governing equations 
are solved iteratively in conical coordinates to yield the complicated wave 
structure of the internal corner and the simple peripheral shock of the 
external corner. The numerical results for the Internal flows compare favor- 
ably with existing experimental data. 

INTRODUCTION 

Existing supersonic aircraft such as the B-1, F-14, F-15., Concorde, and . 
recently proposed designs of advanced hypersonic research aircraft such as 
that shown in figure 1 , possess engine inlets which are composed of planar 
compression or expansion surfaces with swept and unswept leading edges. 

These surfaces form various combinations of Internal and external axial cor- 
ners. Such a corner configuration can generate a rather complicated inter- 
ference flowfield whose prediction is of considerable Interest to the vehicle 
designer, because of the severe pressure gradients and high local heating that 
can occur at the surface. 

The typical internal corner configuration studied in this paper and the 
coordinate system used are shown in figure 2. The flow direction is aligned 
with the x-axls. The vertical wedge is unswept and is always considered a 
compression surface; i.e., 62 > 0. The horizontal or base wedge can have a 
sweep of A and be either a compression or an expansion surface; i.e., 

61 ^ 0. 

The conical wave structure for a typical swept, compression-compression 
configuration is also shown in figure 2. The shock structure consists of a 
planar shock emanating from the leading edge of each wedge, a corner shock 
that joins the two wedge shocks, and two embedded shocks that stretch from 
the body to the triple points. A slip surface exists between each of the 
triple points and the axial corner at which there exists a vortical singu- 
larity. 
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If the base wedge is an expansion surface; i.e., 6 ]^ > 0, then a Prandtl- 
Meyer expansion fan emanates from the leading edge and Intersects the vertical 
wedge shock perpendicularly. As a -result of the Interaction, the wedge shock 
Is bent and shifted ‘slightly towards the axial corner, becoming weakened In 
the process. The expansion fan, after Interaction with the wedge shock, 
becomes concave with respect to the comer. 

A typical external comer configuration is shown in figure 3 and con- 
sists of two Intersecting compression surfaces with angles 5]^ and 62 
sweeps and A£. The flow structure about the external corner Is consid- 

erably simpler than that of the Internal corner. It consists of a single 
continuous shock wave whose strength In rounding the corner transitions from 
that of the d^-wedge shock to that of the 62 -wedge shock. Like the inter- 
nal'-flow problem, there also exists a vortical singularity at the axial cor- 
ner due to the convergence of the cross-flow streamlines (each with different 
values of entropy) . 

Both the Internal and external axial corner flows are conical since 
there is no characteristic length associated with the body. Viscous effects 
are assumed to be minimal, and therefore the flow is governed by the three- 
dimensional steady-flow Euler equations. These equations, under a nonorthog- 
onal coordinate transformation which Introduces conical self -similarity and 
aligns certain Independent variables with the body, are hyperbolic and can be 
Integrated In an Iterative fashion using MacCormack’s finite-difference 
algorithm (ref. 1). The Internal corner problem with Its complicated wave 
structure is solved using the shock-capturing philosophy (ref. 2 ) while the 
external corner problem with its single peripheral shock is solved using both 
the shock-capturing and shock-fitting approaches. 


THEORY 


The governing partial differential equations (continuity, x, y, and z 
momentum) In Cartesian coordinates (see figs. 2 and 3) are first written In 
the dimensionless strong conservation-law form: 
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and where k = (y-1)/2y» Y being the ratio of specific heats. 
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In equation (1), pressure p and density p are made dimensionless with 
respect to free-stream stagnation conditions, and velocity components u, v, 
and w are made dimensionless with respect to the maximxim adiabatic velocity. 
The above system of equations is made complete by the addition of the steady- 
flow energy equation which can be written in the following dimensionless form: 

p = p^l - u^ - v^ - w^) (2) 

It is desirable when solving fluid flow problems to transform the 
Independent variables so that the new coordinates are aligned with the; sur- 
face of the body. This alleviates the numerical difficulties associated with 
satisfying the tangency conditions for an unequally spaced grid. A transfor- 
mation which satisfies the above criterion and also Includes conical self- 
slmllarlty Is 


C = X 

n = y/ (x - z tan A^) 
5 = z/(x + y tan A 2 ) 


(3) 


where A^ and A 2 are depicted In figure 3 for the external corner. For 
the Internal corner, A^ = A and A 2 = 0. 

Under the above transformation equation (1) becomes 


* * * * . . V 

E + F +G_ + H = 0 (4) 
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F = -E pc + F c + G nc + tan 

A AA AA AA 

G = -E 5d - F tan A 2 + G d 

H* = E*[-c(-b +(C/e))tan - d(a -(n/e))tan A 2 ] 

+ F [-d(-l +(l/e))tan A 2 - 5(c/e)tan tan A 2 ] 
+ G [-c(l -(l/e))tan A^^ - n(d/e)tan Aj^ tan A 2 ] 
a = n(l - C tan Aj^)/e_ 
b = 5(1 + n tan A 2 )/e 
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c = C(1 + a- tan A 2 ) 
d = C(1 - b tan A^) 
e, = 1 + n? tan A^^ tan A 2 



Thus equation (4) governs the flow over both the internal and external cor- 
ners and can be solved without any further relations using the shock- ^ 
capturing approach. 

In order to treat the peripheral shock of the external problem as a 
sharp discontinuity, equation (4) is first normalized between the surface of 
each wedge and the shock. This requires a separate coordinate transformation 
for each of the two regions outlined by points 1-2-3-4-1 and 5-6-7 -8-5 in 
figure 4. The necessary equations will be developed for the horizontal 
region only, since the derivation of the analogous equations for the vertical 
region is the same. 


The required transformation is 


1 = (n - n^)/[hg(?,0 “ 


5 = 5 

• I 

and when applied to equation (4) yields 
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The equation of the shock is yg = yg(x,z), and the quantity (hs)^ of 
equation (6) can be expressed as a function of the derivatives (7s) x 

(yg)^ aa follows: 
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The quantity (hg)^ is evaluated numerically and (7s)x computed from 
the ^nklne-Hugonlot equations: 
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and the subscript °° indicates the free-stream condition. 

Following the approach developed by Thomas et al. (ref. 3), it is only 
necessary to know the pressure behind the shock in order to propagate it. The 
remaining flow quantities can be found from equation (8). Shock pressure is 
determined by a finite-difference approximation of the governing partial dif- 
ferential equations .at the shock (refs. 4 and 5).- 

Equatlon (4) , which governs both corner flows (via shock-capturing), and 
equation (6), which governs only the external corner flow (via shock-fitting), 
are hyperbolic with respect to the conical coordinate ?. These equations 
can therefore be solved Iteratively until. E* of equation (4) and of 

equation (6) are zero, indicating the establishment of conical flow. 

The boundary condition at the surface of each wedge requires that .the 
flow be tangent to it. Since an Iterative procedure is employed to .solve the 
governing equations, a simple Euler predictor/modified Euler corrector with, 
one-sided normal derivatives is used at each surface. The following condi- 
tion on the velocity components is imposed after the corrector, step. to. satisfy 
tangency: 
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V = u tan 6 ^ - w tan 6 ^ tan A 2 ; horizontal wedge 
w - u tan 62 + V tan 62 tan A 2 ; vertical wedge 

The axial corner is treated as a multiple-valued point to account for the 
vortical singularity that exists there. Although pressure is continuous, 
density and velocity are discontinuous. 

The computational boundaries for both the internal and external corner 
flows are shown in figure 4. Region I corresponds to the unifom free stream, 
region II to the horizontal wedge flow, region III to the vertical wedge 
flow, and region IV represents the unknown conical flow. Regions 1, II, and 
III are known flow regions and are hyperbolic zones -in the cross-flow planes. 
Region IV is unknown and elliptic in the cross-flow plane. It is the deter- 
mination of this region that is the crux of the present problem. 

RESULTS 


Numerical results for the internal corner flow problem are shown in 
figures 5-10. These results were computed by Kutler (ref. 6 ) and Shankar and 
Anderson (ref. 7). 

Some of the most recently published experimental data obtained for the 
corner-flow problem are by West and Korkegi (ref . 8 ) . They tested an equal 
wedge angle ( 6 ^ ~ &2 ~ 9.49") configuration in Mach 2.98 flow over a 
Reynolds number range from 0.4x10® to 60x10®, which Included laminar, tran- 
sitional, and turbulent boundary layers. A numerical solution for this same 
case was obtained, and the shock wave and slip surface structure are compared 
with the high Reynolds number experiment in figure 5. 

The inviscid embedded shock is slightly concave when viewed from the 
origin falling inside the location' of the corresponding experimental shock. 
The corner shock, which is slightly convexed when viewed from the origin, 
also falls inside the experimental shock. The positions of the experimental 
. and numerical wedge shocks agree exactly. It appears, therefore, that the 
displacement effects of the boundary layer in the region bounded by the cor-r 
ner and embedded shocks result in an effective thickening of the body, and 
this forces the shock structure outward. 

The location of the slip surfaces for this case can be found from plots 
of density and is shown as the thin double line in figure 5 stretching from 
the triple point to the origin. The slip surface is slightly curved and 
asymptotically approaches the bisector near the origin. The experimental 
shear layer is also curved but appears to merge before the origin is reached. 
Since the positions of the numerical and experimental triple points are dif- 
ferent, the comparison between the inviscid slip surface and vihcous shear 
layer, which originate at the triple points, is unfair. But, qualitatively, 
their basic -shapes are the same. - 

A comparison of the numerical and experimental (turbulent boundary 
layer) surface pressures is shown in figure 6 . The first pressure rise in 
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the experimental data. (decreasing) indicates the onset of separation. This 
is followed by a reduced gradient region that indicates separatlpn and again 
a rapid pressure rise that indicates reattachment. The pressure between the 
reattachment point and the. origin is greater than that of. the inviscid " 

result. This higher pressure indicates an apparent thickening of the body in 
this region due tp -boundary-layer displacement effects. 

Nangla (ref. 9) performed an experimental study on the wave interactions 
in supersonic intakes and obtained some rather interesting data. Figures 7-9 
are numerical solutions which compare with Nangla 's results. Figure 7 is a 
pressure contour plot for an unswept configuration with 6 ]^ = -5® and 
62 =7.5“ for a Mach number of 3. It is interesting to note that the wave., 
structure does not exhibit any corner shock or slip surfaces. This appears 
to be correct because changes through the expansion fan occur Isentroplcally . 
However, an auxiliary compression wave is formed in the corner region. The ' 
wedge expansion fan, after the interaction with the wedge shock, turns away 
from the axial corner and is concave when viewed from the origin, whereas the 
wedge shock turns toward the axial corner and appears to be slightly convex ^ . 

The inviscid numerical wave structure compares very well with the exper- 
imental results of Nangla. The auxiliary compression fan in the corner 
region, as predicted by the numerical solution, is not observed in the exper-' 
Imental results, however. 

A pressure contour plot for a A = 30“, 6 ]^ f 5“ , and 62 = 7.5“ .cohfig- 
j uratlon at. a Mach number of 3 is shown in figure 8 . The numerical solution 
is again compared with the experimental data of Nangla. The computed surface- 
pressure distribution for both the horizontal and vertical wedges is compared 
with Nangla 's experimental measurements in figure 9. ... .V" , 

The final case considered on the Internal corner problem consists of a ’ ^ 
A =:30“; 61 = -5“, and 62 = 7.5“ configuration at a Mach number of 3... The 
• pressure contour plot of the numerical solution is shown in figure 10. The. 
wedge expansion fan, after interacting with the wedge shpckj. turns away from 
the axial corner and appears to be concave when viewed from the origin. The 
intersecting wedge shock is deflected towards the axial corner under the 
influence of the wedge expansion fan.. An auxiliary weak shock is formed in 
the corner region which merges with the weak embedded shock and forms a 
strong. shock near the horizontal wedge surface. The shock impingement on 
this surface is nearer the axial corner than the wedge shock. The strength 
of the impinging shock appears to increase with increasing sweep angle.. 
Furthermore, the Interference region decreases as the sweep single increases.. 

Numerical solutions for two external corner configurations were generated 
and are shown in figures 11-14. Results are presented for both the shock- 
capturing and., shock-fitting techniques. Figure 11 shows the shock, shape and* 
cross^flow sonic line for an unswept equal-angled ( 6 j^ = 62 = 10 “) .configura- 
tion at a Mach number of 3. The surface-pressure distribution for, this case 
compared to linear theory is shown in figure 12 . 

The shock and cross-flow sonic-line locations for a swept (A^ = A 2 = 

30“)., equal wedge angle (61 = 62 = 10“) configuration at Mach nximber 3 is 


650 



shown In figure 13. The shock-capturing results are almost identical to the 
shock-fitting solution. The surface pressure distribution for this configura- 
tion is shown figure 14, and again both the shock-capturing and shock-fitting 
results agree . 

All of the numerical results presented here were obtained oh serial 
machines but employed a numerical procedure which is particularly well suited 
for thfB parallel processing philosophy. Should the computing power of an 
ILLIAC IV be required for a more complicated problem, such as a multiple- 
comer configuration, some of the coding techniques used for the solutions 
presented here could easily be applied. 
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Fig. 1 Hypersonic airbreathing aircraft 


y 








■ Fig. 5 Comparison of numerical and experimental shock patterns; M « 2.98, 
6 l = 62 == 9.49“. 



Fig. 6 Comparison of numerical and experimental surface-pressure distribu- 
tion; M = 2.98, 61 - 62 = 9.49, A = 0“. 


654 



Fig. 7 Comparison of numerical and experimental wave structure; M =-3, 
6i = -5“, «2 = 7.5% A = 0“. 



Fig .*' 8 Comparison of numerical and experimental shock structure; M = 3, 
6 i = 5°, 62 = 7.5“, A = 30®. 
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Fig. 9 Comparison of numerical and experimental surface-pressure distribu- 
tion; M = 3, 6 i =5% 62 = 7.5", A = 30". 



Fig. 10 Pressure contour plot of computational plane; M = 3 , = -5", 

62 = 7.5", A = 30". 
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Fig. 11 Shock structure for S 3 nnmetrlc external corner; M = 3, 5l = (S 2 “ 10 
Ai = A2 = 0®. 
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Fig I 13 Shock structure for swept external corner; M = 3, 6 i “ ^2 “10 
Ai = A2 = 30®. 
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Fig. 14 . Surface-pressure distribution for swept external corner; M 
6 i = 69 = 10°. Ai = A 2 = 30*. . . 



